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Abstract . A simple heuristic argument is presented which suggests that 
the large rigidity independent mean free paths observed for low rigidity 
solar and galactic cosmic rays can be understood in terms of weak turbu- 
lence diffusion theory if one assumes that the interplanetary magnetic 
turbulence consists of a combination of Alfven waves propagating with 
constant field magnitude and a small (</‘5-10K) admixture of compressive 
fluctuations. 


1 . Introduction 

In a recent series of papers Zwickel and Webber [1977, 1978] and Fisk 
[1979] have discussed what appear to be serious discrepancies between the 
theory of cosmic ray propagation and observations of low rigidity solar 
particles. In the first place, the mean free path, X, for pitch angle 
scattering computed from weak turbulence theory, and the value inferred 
from solar particle events differ by more than an order of magnitude. In 
addition, at low rigidities X is observed to be approximately constant with 
decreasing rigidity, P, below a few GV, while from theory one expects X to 
decrease. The essence of the problem is apparent in Figure 1, taken from 
Fisk [1979]. The quasilinear resonant scattering theory [Jokipii, 1971] 
yields the solid curve denoted "theoretical upper limit." Nonlinear 
corrections to quasilinear theory [Volk, 1975; Goldstein, 1976; Jones et 


al., 1978] worsen the discrepancy because of the enchanced scattering at 
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pitch angles, a, near 90® (a is the particle's pitch 2 uigle defined with 
respect to the mean field, V- 

The purpose of this letter is to suggest a simple physical argument 
which yields both the correct rigidity dependence and magnitude of x. In 
the following section we briefly review those aspects of current theory 
leading to the prediction of very small mean free paths at low rigidities. 
We then argue that observations of interplanetary MHD waves indicate that 
the character of the magnetic turbulence has not been sufficiently well 
represented by current analyses. These observations suggest a simple 
quantitative resolution to the discrepancy. 

2. Origin of the problem 

The theoretical result that X decreases with decreasing P, arises 
primarily from the very strong cyclotron resonance between charged 
particles and parallel propagating Alfven waves (parallel being defined 
v.ith respect to The inclusion of compressive MHD modes to the 
turbulence only enhances the wave-particle scattering because of the 
resultant strong Landau resonance [Goldstein et al., 1975; Lee and Volk, 
1975]. Only in the case of turbulence composed of wavevectors, jt, 
propagating predominately at angles nearly 60® to the direction of will 
pitch angle scattering be inefficient in backscattering particles. This 
possibility has been developed in some detail by Morfill et al. [1976], 
but as far as is presently known the interplanetary turbulence consists 
primarily of Alfven waves propagating significantly closer to the direction 
of than expected from WKB analysis [Denskat and Burlaga, 1977]. The 
observations are by no means exhaustive, and the suggestions of Morfill et 
al. [1976] may well prove Important. In addition, the rigidity dependence 
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cannot easily be accounted for in that picture. For that reason, we feel 
that the basic explanation of the observed mean free paths lies elesewhere. 

There is at least one area where the weak turbulence theory has proved 
successful. That is in the apparently excellent agreement between 
numerical simulations of the pitch angle scattering coefficient, 
(m=coso), conducted by Kaiser et al. [1978] and the predictions of the 
nonlinear theories [Goldstein, 1976; Jones et al., 1978]. In those 
comparisons the theories were carefully tailored to the simulations in that 
the power spectrum of magnetic turbulence was nearly identical in both 
treatments. In addition, the simulations were designed to minimize the 
effects of transients, and therefore modeled a purely steady state and 
diffusi/e situation. Similarly, because the adiabatic approximation is 
made in all the nonlinear theories, they also describe diffusive propaga- 
tion. In a recent series of papers, Klimas et al. [1976a, b, 1977] have 
argued that the adiabatic approximation cannot be made, and that the 
Hfinadiabatic, or diabatic, behavior of the initial transients plays a 
crucial role in the long time relaxation of the particle distribution 
function toward isotropy. Thus a problem may exist in relating the 
simulations and nonlinear theories to interplanetary observation?., ..’here a 
steady state cannot be assumed and where the role of diabatic effects may 
be crucial. Unfortunately, the analysis of Klimas et al. was confined to 
the use of quasilinear, helical, particle orbits, which they showed to be 
inadequate. Because their theory has not been extended to include perturb- 
ed orbit effects the consequences of their arguments are presently unknown. 
With this reservation in mind, we will proceed with our discussion, taking 
the agreement between simulations and nonlinear theories as an indication 
that the theories do approximately describe charged particle propagation, 
and that the language and ideas of diffusion theory are appropriate for the 
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interplanetary problem. 

When one attempts co use the correlation function employed in the 
simulations as an approximation to the interplanetary turbulence, the 
relative ease with which particles can propagate through a=90® (which 
resulted in the close agreement between theory and numerical experiment) 
produces very small values of X. The relationship between and the 
spatial diffusion coefficient, <||, from which one finds X is given by 
' Hasselmann and Wibberenz, 1968; Earl, 1974] 

1 1 (1 - 

•'ll = _ /dw (1) 

M 2 » (2D^) 

(The dimensicnless notation follows Goldstein [1976 and 1977] where times 
are measured in gyroperiods and lengths in Larmor radii.) In the nonlinear 
theories, strong nonlinear perturbation of the particle trajectories 
produce nearly isotropic scattering, so that one could have just as well 
usea the expression [e.g. Jokipii, 1971; Goldstein, 1977] 



The physics of the scattering near p=0 has been analyzed by VoiK [1975] 
-ind Jones et al. [1978], who showed it is due almost entirely to particles 
mirroring in tne random fields. The mirroring arises even in the case of 
turbulence consisting of linearly or circularly polarized Alfven waves (or 
tne eu'jivalent magnetostatic slab model) because of a second order change 
1 .. field magnitude. This interpretation is also consistent with the 
scaling laws found for D^(u=0) when e=(X^/rg)>1 (low rigidities); namely, 

y 

cli^(u=0)«n with v»3, approximately independent of e. (X^ is the correla- 
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tion length of the magnetic turbulence, and is the particle's Larmor 
radius.) Goldstein et al. [‘'975] have pointed out that for eD^(u=0) is 
independent of rigidity when mirroring predominates. One should note that 
for wave fields having the property that the total magnitude of fi. is a 
constant, the theories of Jones et al. [1978] and Volk [1975] would predict 
virtually no scattering through usO . Such wave fields are, for example, 
characteristic of Alfven waves observed in the solar wind [Burlaga and 
Turner, 1976]. 

Goldstein [1976] gave a different scaling at m= 0 which resulted from the 
momentum dependent terms absent in either of the other nonlinear theories. 
In that case, the scaling became €D^(M=0)«n'^e for t>^ . Thus, there are 
additional physical effects tending to backscatter particles. However, at 
the low rigidities at which these momentum dependent terms become impor- 
tant, e is so large that the validity of the weak tubulence theory itself 
is questionable [Goldstein, 1977]. 

3. The mean free path at low rigidities 

Our proposed explanation for the observed behavior of low rigidity 
cosmic rays is based on our contention that the models of interplanetary 
turbulence thus far employed in simulations and theory are inaccurate in at 
least one very important respect. 

Observations of interplanetary MHD waves by Burlaga and Turner [1976] 
have shown that the root mean square value of the fluctuation in |£| is 
about 0.06 IB.I (or n=0.06). In such a wave field, the amount of second 
order mirroring present and capable of scattering particles through p=0 is 
very limited. Thus we can conclude that both the extant numerical 
simulations and the theoretical analyses greatly overestimate the amount of 
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backscattering in the solar wind because all these models use correlation 
functions which dc not preserve field magnitude in second order. This 
brings us to our basic hypothesis, whose implications we develop below; 
the xaterulanetarv Alfvenlc wave tu rbulence conserves |B| and that residual 


Obviously, in the vicinity 
of fast and slow stream interaction regions where |Ii| changes, this 
requires modification. There the value of n for compressive modes probably 
exceeds 0.1 and particles encountering tnese inter'action regions should be 
easily scattered [cf. Barouch and Burlaga, 1976] which seems consistent 
with recent observations of Gold and Roelof [1979]. 

If „ne Alivenic tui buienoe conserves |£.| , tnen in tne absence of any 
other wave modes there should be greatly inhibited propagation through p=0. 
(Several years ago Sari, ^.n an unpublished preprint discussed a similar 
suggestion.) In fact, D^(ii=0) should be approximately given by the 
quasilinear formulas, and the strong, nearly isotropic scattering seen in 
the current simulations should be absent. If one now takes into account 
Derskat and Burlaga' o [1977] observation that the angle, 0, between and 
is in gancral small but nonzero, then from the well known properties of 
quasilinear theory we conclude tnat with X+® i''or interplanetary 

Alivenic turbuience [sec Volk, 1975 and Jones et al., 1973]. Recently, 
Kliraas et al. [1976] nave claimed '•hat even f?'r 0=0 the quasilinear 
(adiabatic) theory does not allow particle propagation through u=0. 
Furthermore, when the adiabatic approximation is not made, Klimas et al. 
[KTSc.b, ■’9771 showed that particle propagation through u=0 is greatly 
i.-.hib.ited , with a iogarithmic relaxation ol the particle distribution 
t diction toward isotropy 

Thus we are led to the following ansatz: In the presence of magnetic 

turbulence consisting only of Alfvenic fluctuations conserving field 
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magnitude, particle propagation through u=0 is essentially zero. We now 
proceed to explore the consequences of this ansatz, leaving for future work 
the task of verifying it through more formal analysis. We go on to 
consider the residual effect of the 6$ change in |fl.| in determining both 
the magnitude and rigidity dependence of X. 

In the magnetostatic approximation, fluctuations in |fi.| give rise to a 
Landau resonance of the form D^(p=0)«6(u) [Fisk et al., 1974; Goldstein et 
al., 1975; for a generalization to finite frequencies see Lee and Volk, 
1975]. Volk [1975] has argued that the effect of nonlinear perturbations 
is to spread the Landau resonance over a finite region of y-space while 
preserving its integrated weight. In that case, because D^( compressive 
modes) is maximum at y=0, and non-zero within the range |y|<±0.3, equation 
(2) can be used to calculate '^j|* This has been done by Klimas and Sandri 
[1973]. Assuming n=0.1, we find that at low rigidities X is independent of 
rigidity with a value of approximately 0.3 AU (cf. Fig. 3 of Klimas and 
Sandri [1973]). Note that Klimas and Sandri's calculation included the 
effects of some scattering by Alfven waves in that they used a correlation 
function typical of isotropic turbulence. However, the Landau resonance i' 
so efficient in scattering particles near y=0 that the cyclotron resonances 
can be completely ignored. This statement can be strengthened further by 
noting that because the observed Alfven waves do not propagate precisely 
along fi, the efficiency of the cyclotron resonance in even bringing 
particles close to y=0 is further reduced. 

At high rigidities, e<<1, resonant scattering will again become 
important and the standard results should obtain using the full rms value 
of n“0.3 with X^»2«10‘® cm [Fisk and Sari, 1973]. The composite theoreti- 
cal curve is shown in Figure 1 . 
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4. Discussion and conclusions 

The essential ideas of this letter can now be summarized. Alfvenic or 
slab turbulence conserving l£| will not appreciably backscatter low 
rigidity particles. The residual 5-10$ fluctuations in |£.| give rise to a 
value for X which nas both the correct magnitude and rigidity dependence. 

These ideas are amenable to test by two techniques. Numerical simulations 
can be constructed in whch the turbulence conserves |£,| [Owens, provate 
ccmnunication] . Furthermore, the nonlinear theories of Volk [1975], 
Goldstein [1976] and Jones et al. [1978] can be modified to use correlation 
functions subject to that same constraint. 

A ocnowlc-lgements. it is a pleasure to thank Drs. A. J. Klimas, F. C. 

Jones and L. F. Burlaga for many stimulating and illuminating discussions, 
ana Dr. J. R. Jokipii for stimulating my renewed interest in this problem. 
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Figure Caption 


Representative values of X|| deduced from solar cosmic ray observations. 
The solid line labeled ”Theoretical upper limit” is from weak turbulence 
calculations using equation (1) for iC|| and correlation functions previously 
thought to model interplanetary magnetic turbulence. The dashed curve is our 


estimate of the residual scattering assuming that the Alfvenic turbulence 
conserves |£,1 and all scattering results from the residual compressive 


fluctuations (from Zwickel and Webber [19773 and Fisk [1979]). 



McKibben et al. 



Rigidity (GV) 





BIBLIOGRAPHIC DATA SHEET 


1. Report No. 2. Government Accession No. 

TM 80290 

3. Recipiint’s Catalog No. 

4. Title and Subtitle 

The Mean Free Path of Low Rigidity Cosmic Rays 

5. Report Date 
May 1979 

6. Performing Organization Code 

7. Author(s) 

Melvyn L. Goldstein 

8. Performing Organization Report No. 

9. Performing Organization Name and Address 
NAS A/ GSFC 

Laboratory for Extraterrestrial Physics 
Interplanetary Physics Branch, Code 692 
Greenbelt, MD 20771 

10. Work Unit No. 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

1 2. Sponsoring Agency Name and Address 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 


A simple heuristic argument is presented which suggests that the large 
rigidity independent mean free paths observed for low rigidity solar and 
galactic cosmic rays can be understood in terms of weak turbulence diffusion 
theory if one assumes that the interplanetary magnetic turbulence consists 
of a combination of Alfven waves propagat ng with constant field magnitude 
and a small (^^5-10%) admixture of compressive fluctuations. 



19. Security Classif. (of thi% report) 20. Security Classif. (ot this page) 21 No. of Pages 22. Puce* 


U r 13 


For sale by the National Technical Intormalmn Service, Springfield Virginia ?7tb1 


GSFC 25 44 (10/77) 






















